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We calculate the surface critical exponents of the ordinary transition occuring in semi-infinite, 
quenched dilute Ising-like systems. This is done by applying the field theoretic approach directly in 
d = 3 dimensions up to the two-loop approximation, as well as in d = 4 — e dimensions. At d = 4 — e 
we extend, up to the next-to-leading order, the previous first-order results of the ^/e expansion 
by Ohno and Okabe [Phys. Rev. B 46, 5917 (1992)]. In both cases the numerical estimates for 
surface exponents are computed using Pade approximants extrapolating the perturbation theory 
expansions. The obtained results indicate that the critical behavior of semi-infinite systems with 
quenched bulk disorder is characterized by the new set of surface critical exponents. 

I. INTRODUCTION 

The critical behavior of quenched random systems undergoing continuous phase transitions is of great interest. 
It is well known that the critical behavior of ideal pure bulk substances may be changed by introducing disorder. 
A prominent result in the theory of quenched disordered systems is the Harris criterion [Q which states that the 
presence of disorder is relevant for those pure systems which have a positive specific heat exponent a. Thus, in the 
class of 0{N) symmetric iV-vector models in d space dimensions the Ising model is the one of primary interest, having 
a{d) > 0. The marginal value, a = 0, corresponds to d = 2. In this case the small amount of disorder produces a 
marginal perturbation. This results in logarithmic corrections to the power-law singularities without modification of 
the critical exponents with respect to the pure system (for reviews see |^J^). 

In three-dimensional disordered Ising systems a new set of modified critical exponents appears. This was confirmed 
by renormalization group (RG) calculations 0-0], experiments Ppl, and Monte-Carlo simulations [^0. After 
pioneering investigations using the Wilson's RG and e expansion |ffl4 l2| , |l3| , scaling- field method |Q and the massive 
field theory in three dimensions )p"5|-p^, there has been a great renewed interest to the subject in the last few years 
p9|-p3[. An interesting relation between the random Ising model and "iV-colored" tethered membranes has been 
given in Various aspects of physics of quenched disordered systems are the subject of recent extensive MC 

investigations |2|-|2^. 

The presence of boundaries, which are inevitable in real systems, leads to the new, surface physics. General reviews 
on surface critical phenomena are given in |8|-|30|. It is now well known that there are several surface universality 



classes defining the critical behavior in the vicinity of boundaries, at the temperatures close to the bulk critical point. 
Each surface universality class is defined both by the bulk universality class and specific properties of a given boundary. 

What happens with the surface critical behavior, when the quenched disorder is introduced? First, the defects may 
be localized only at the boundary. In pl[ , the relevance- irrelevance criteria of the Harris type concerning the weak 
surface disorder have been worked out. In the case of the ordinary transition (which corresponds to the free boundary 
conditions), it has been demonstarated that for d > 2 the weak surface randomness is an irrelevant perturbation. The 
same conclusion follows from rigorous arguments ]^^ and Monte Carlo simulations js^. There is a strong evidence 
that a small amount of quenched surface disorder is irrelevant at the special transition as well |^,^,|^ . 

Second, the random quenched disorder may be introduced in the bulk, say, of a semi-infinite system bounded by a 
plane surface. This will, in general, shift the critical temperature of the bulk phase transition, and drive the system 
to another, "random" fixed point. According to the Harris criterion, the Ising systems are of main concern here. The 
change of the universality class of the bulk will affect the critical behavior of the bounding surface. Thus one will 
expect the new surface critical exponents to appear. And, in view of the aforementioned properties of boundaries, 
these new surface exponents should emerge irrespective of the presence of the disorder just at the boundary itself. 

This is a situation typical for d > 2. The problem of a RG calculation of "disorder-induced" surface critical 
exponents was first addressed by Ohno and Okabe in 1992 by using the e expansion about the upper critical 
dimension d = A. In the marginal case d ~ 2, the boundary critical behavior of random systems has been studied 
very recently in a series of papers ||3^-|3^ . 

In the present paper, we present an attempt to go beyond the approaches and approximations employed in p^ ] 
in calculating the surface critical exponents of random semi-infinite Ising-like models. Our main caclulations are 
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performed directly at d = 3 using the massive field theory approach [^9|-^2|j3J] in the two-loop approximation. 
Moreover, we extend, up to the next-to leading order of the -y/g expansion, the previous first-order results of Ohno 
and Okabe |^^. In both cases, numerical estimates for surface critical exponents of the ordinary transition are 
evaluated using Pade approximants improving the direct perturbation theory expansions. The obtained values of the 
surface critical exponents are consistent with the results by Ohno and Okabe jSSt and confirm that the semi-infinite 
systems with quenched bulk disorder are characterized by a new set of the surface critical exponents. 



II. MODEL 



The description of the surface critical behavior at the ordinary transition can be formulated in terms of the effective 
Landau-Ginzburg- Wilson Hamiltonian 



dzld 



(2.1) 



The d-dimensional spatial integration is extended over a half-space M'^. = {x={r, z) G M"^ \ r G ]R'^~^,z > 0} bounded 
by a plane surface at z = 0. Here f = ip{x) is a continuous scalar field corresponding to the one-component order 
parameter of an original Ising system. The surface is considered to be free, and hence the field if{x) satisfies Dirichlet 
boundary conditions at z = 0, ^(r,0) = |4^j29[ |. 

One of the possibilities to introduce disorder into the model is to assume that the parameter tq incorporates local 
random temperature fluctuations 6t{x) via tq = mQ + dT{x). Here ttiq is a "bare mass" representing linear temperature 
deviations from the mean-field ctitical temperature, and the random variable St(x) has the properties ((5T(x))conf = 
and ((5T(x)(5T(a;'))conf — A(5(a: — x') with A > 0. Angular brackets with the subscript con/ denote configurational 
averaging over quenched disorder which should be implemented on the level of the free energy plJi^ ] . 

Ohno and Okabe | |35[ | analyzed the above random model by direct averaging over disorder using the method originally 
introduced by Lubensky |^ . Alternatively, the configurational averaging of the free energy can be performed using 
the replica trick 



F = -T hm -((Z")eo„f-l) 



(2.2) 



where Z is the partition function of a configuration given by the Boltzman weig ht e-HI-^-l, as it was first done in the 
RG calculations by Grinstcin and Luther |Q. We shall use this last possibility to treat randomness. Employing the 
replica trick leads to the effective LGW Hamiltonian with cubic anisotropy defined in the semi-infinite space, 



dz\d'^-^i 



(2.3) 



where (p{x) is an n- vector field with the components (pi{x), i = 1, • • • , n. This last effective Hamiltonian is appropriate 
for the description of the ordinary trans itio n of random systems in the replica limit n — > 0. The 0{n) symmetric term 
arises due to the random averaging in (2.2) via cumulant expansion and thus its coupling constant uq cx —A < 0. 



III. CORRELATION FUNCTIONS AND THEIR RENORMALIZATIONS 



The model defined in (2.2) is translationally invariant in directions parallel to the boundary. Thus it is often useful 
to perform Fourier transformations in d — 1 dimensional subspace with respect to "parallel" coordinates r. We shall 
denote the associated parallel momenta as p. In the perpendicular direction the coordinate z is retained. Sometimes, 
in the perturbative calculations it is advantageous to work in the complete coordinate representation, without any 
transformations to the momentum space. This is another difference with respect to the approach of Ohno and Okabe 
who employ the full momentum-space representation in all d directions as it was originally done in the earliest 
field-theoretical work on semi- infinite systems BSMffl. 



The fundamental two-point correlation function of the free theory corresponding to (2^) is given by the Dirichlet 
propagator: 

{(pi{r,z)ipj{0,z'))o = Gr,ir;z,z')Sij . (3.1) 
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In the familiar pz representation the Dirichlet propagator reads 

1 
2ko 



Gb(p; z, z') = [e-"°l"-"'l - e-^"^^'+''A , (3.2) 



where the standard notation is used, Kq — \JjP' + rn^ . 

The Dirichlet propagator vanishes identically when at least one of its z coordinates is zero. Consequently, all the 
correlation functions involving at least one field at the surface vanish. This property holds for both the free and 
renormalized theories pOt . Owing to this property, it is not a straightforward matter to analyze the surface critical 
behavior at the ordinary transition. 

In fact, the critical surface singularities at the ordinary transition can be extracted by stud ying t he correlation 
functions with insertions of (inner) normal derivatives of the fields at the boundary, dn4>{r) ]4^ ^ ^ , ^^ . Actually, in 
order to obtain the characteristic exponent ri°i^^'^ of surface correlations, it is sufficient to consider a correlation function 
of two normal derivatives of boundary fields 

S2{P) = ^{p,z)\.=o ^ ^i~P,^')Uo) ^ (3-3) 

where the fields ^p(j>, z) are the Fourier transforms of the fields ip{r, z) in rf — 1 dimensional parallel subspace. Q2{p) is 
a parallel Fourier transform of the corresponding two-point function Q2 {r) in direct space. At the critical point Q2 (p) 

behaves as p~^~^^ii . It reproduces the leading critical behavior of a two-point function G2{p) = {(p{p, z)(p{—p, z )) 
in the vicinity of the boundary plane. The surface critical exponent 7]°^^"^ is given by the scaling dimension of the 
boundary operator dn^p{r). 

In the presence of randomness, the exponent r]°^'^ differs from its counterpart in ordered semi-infinite system. The 
other surface exponents of the ordinary transition can be determined through the scaling laws [p9| . 

In the present formulation of the problem, the renormalization process for the random system is essentially the 
same as in the "pure" case [ p9|j3^ . One introduces the renormalized bulk field and its normal derivative at the surface 
through 

iPb{x)^Z^^ tf{x) and (9„(p(r))^=Z|| ^5„(^(r), (3.4) 

and renormalized correlation functions involving N bulk fields and M normal derivatives, 

g^j^^''\p-m,u,v) = Z-'^Z-^g'^^'^'\p-mo,uo,vo) (3.5) 

for {N,M) (0,2). In the special case of only two surface operators {N,M) = (0,2), an additional, additive 
renormalization is required, so that 

g^°-'^ ip) = z-^' _ ^(0,2) (p^o)] (3 g) 

Everywhere the renormalizations of the mass and coupling constants are implicit. These are the standard ones of 
the massive infinite- volume theory. The surface renormalization factor Z\\ (u, v) can be conveniently obtained from 
the consideration of the boundary two-point function which is the object with the simplest Feynman graph 

expansion (see below): 

Z||=-limI^^C?(°'^)(p). (3.7) 
p^o p c)p 

A standard RG argument involving an inhomogeneous Callan-Symanzik equation yields the anomalous dimension of 
the operator 9„(p(r) 

d 

li^^ll I FP (3.8) 

ainZ||(u,w) ainZ||(M,u) 
= Pu[u,v) h Pv[U,V) 



du dv 



FP 



"FP" indicates here that the above value should be calculated at the infrared-stable random fixed point of the 
underlying bulk theory, {u,v) = {u*,v*). The surface critical exponent rj°^'^ at the ordinary transition is then given 

by 



= 2 + 7711. (3.9) 
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IV. TWO-LOOP APPROXIMATION 



The Feynman diagram expansion of the unrenormahzed correlation function Q^^'^'> is, to the two-loop order, 

g^°''^\p; mo, uo,vq) = -kq + 




(4.1) 

Full internal lines represent here the free Dirichlet propagators ( [3.2[ ) and dashed external lines give the factors e"**"^' 
when attached to the internal point with the coordinate z^. Let us enumerate diagrams in the above sequence 1,2,3,4. 
In the present theory described by the effective Hamiltonian (|2.3|), these graphs have their corresponding weights 



(0) 



^(0) 
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with 4"^ = - 



with 4") = ^ 



■Uo + Vo, 



,(0) 



and ^ with t(") = t(°) - ^^W^ 



(4.2a) 
(4.2b) 
(4.2c) 



The typical bulk short-distance singularities, which are present in the graphic expansion (4.1), are subtracted after 
performing the mass renormalization 



,_d_ 



(4.3) 



Here the full lines with sighns "-" denote the free bulk propagators. They are associated with the first term in the 
Dirichlet propagator (3_^), which is the usual bulk massive propagator in the pz representation. As a result of the 
mass renormalization one obtains 



1 



+- 



2k 



vn? d 
2^9fc2 




(4.4) 

In this expansion the "surface" divergences still remain present, due to the one-loop self-energy insertions. These 
are represented by the closed lines with the index "-I-" which means only the second, mirror term of the fre e pro pagator 
( |3.2|) . Moreover, we have explicitly indexed here, with the label "D", the original Dirichlet lines (cf. ( [4.l| )). Bulk 
propagators are denoted with "-" sighns. The last two graphs in the second line represent just the usual bulk 
subtractions. 

The boundary singularity in the one-loop diagram of the two-poi nt su rface correlation function has to be removed 
through the additive renormalization [zero-momentum subtraction ( |3.6| )]. This, however, does not influence the cal- 
culation of the renormalization factor Z\\ which involves a momentum differentiation (see (3.7)). Surface divergences, 
present in each of the two last bubble graphs, mutually cancel. Actually, the whole combination in the third line of 
(4.4) vanishes identically as in [p4|. Hence, applying the rule (3.7), we obtain 



Performing the integration of (E^) by analogy with pj] we derive the result 



Z||(uo,i;o) = l + V+4°^^, (4.6) 
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where the constant C stems from the two-loop contribution into (4.5), 

107 7 4 

C ~ — - - In - - 0.094299 ~ -0.105063 . (4.7) 
162 3 3 

Here the renormalization factor Z|| is expressed as a second-order series expansion in powers of bare dimensionless 
parameters uq — uo/(87rm) and % — vo/(87rm). The corresponding weighting factors t^"'' and 1^"^ are obtained by re- 
placements (uq, vq) — > (uo, vo) in the original combinations tf^ and 4°^ from (4^). As it is usual in SM^jerrenormaHzable 
theories, the renormalization factor expressed in terms of unrenormalized coupling constants is finite. 

As a next step, the vertex renormalizations should be carried out. To the present accuracy, they are given by 

n + 8 _ _\ /.ON 
uo^u[l-\ —u + vj, (4.8a) 

vo^v(l+ ''^v + 2uj . (4.8b) 

Again, the vertex renormalization at d = 3 is a finite reparametrization. All relevant singularities have been removed 
already after the mass renormalization and taking into account the special bubble-graph combinations emerging in 
the theory with Dirichlet propagators. The result is a modified series expansion 

, , n+2 _ V n+2 / n+8\ 
3\ _2 , „ . "-+8 



C+-j v' + 2[C+—j uv. (4.9) 

Combining the renormalization factor Z|| (m, v) together with the one-loop pieces of the beta functions 

l3y,{u,v) ^ -u(^l - ^^^^ u- v^, (4.10a) 

P^iu,v) = -v(l~^v~2uy (4.10b) 

through (^) yields the desired series expansion for ryu . 

In terms of renormalized coupling constants u and v, normalized in a standard fashion so that u=^^^^^u and v=^v, 
we obtain finally 

n+2 V 

^ll("'") = ~2(;^"-6 ^'-^'^ 



(n+2 , o 8^, , o 16 
(n+8)2 ^ ' % ^ ' n+8 



where C(n) is a function of the replica number n, defined as 



n + 1 4 

C{n) = C+^. (4.12) 



In fact, the last expression (4.11) for t^h provides a result for the cubic anisotropic model given by the effective 
Hamiltonian (2.3) with general number n of order-parameter components. In the case of infinite space, this last model 
attracted much attention very recently (see e.g. |2^,^,Q and references therein). 

In the following we restrict our discussion to the case of random Ising system by taking the replica limit n ^ 0. 
Hence, we obtain the next two-loop expansion for the surface critical exponent rj?/''^ 



^f'' = 2 - ^ - ^ - ^C(0)^2 - ^C(l)^;^ - 2C(0)uz;. (4.13) 

Through the scaling relations we get access to the other surface critical exponents. For convenience further below we 
suppress the superscript ord at the surface critical exponents. 



5 



V. SURFACE CRITICAL EXPONENTS 



As it is well-known, the knowledge of one certain exponent at the ordinary transition allows one to define the 
complete set of other surface critical exponents through the scaling relations pOt - For convenience we quote them 
here: 



VI- = 



V + V\\ 




2 




-(d-2- 


^V\\), 


= 1^(1 - V\\ 










), 


A dH 


- 2 - r? 


/3i d- 


2 + 7711' 


Ai 


d - 7711 


Pi d 


-2 + 7?|| 



The exponent ri± characterizes the critical-point correlations perpendicular to the surface, (3i describes the decay of 
the surface magnetization on approaching the critical temperature, 711 is the (local) surface susceptibility exponent, 
71 is the layer susceptibility exponent, Ai is the surface magnetic shift exponent, and dn and give relations between 
the surface magnetization and the surface and bulk external magnetic fields, respectively, along the critical isotherm. 
The values ly, rj, and A — v{d + 2 — 77)72 are the standart bulk exponents. 

In order to o btain individual RG expansions for each surface exponent, we use the above scaling laws (with d = 3) 
combining Eq. ( 4.13 ) with the — > limits of the two-loop series expansions for bulk exponents v and 77 16|T5|j5l|| 



V 



_ 1 
~ 2 

1 
324 
_ _8_ 
~ 27 



1 



11 



6 2(n-h8) 



^ ; (2777 - 38) w - "^1"^^^^^ {27n ~ 38)u^ 



(ri + i 



V 

27 



2uv 



in + 2) 



3{n + 8) (n + 8)2 

For each of surface critical exponents we obtain a double series expansion in powers of u and v of the form 



(5.1) 
(5.2) 



(5.3) 



truncated at the second order. Power series expansions of this kind are known to be generally divergent due to a 
nearly factorial growth of expansion coefficients at large orders of perturbation theory |52|-|55t]. Hence, the numerical 
evaluation of the exponents represented by such series expansions requires additional " resummation" procedures. 

The simplest way to obtain meaningful and rather accurate numerical estimates is to construct a table of rational 
approximants in two variables from the original series expansions. This should work already well when the series 
behave in lowest orders " in a convergent fashion" . Apparently divergent ones require more sophisticated summation 
procedures. 

The results of our Fade analysis are represented in Table I. We evaluate the exponents at the standard two-loop 
random fixed point 03] 



-0.60509, 



2.39631. 



(5.4) 



To give an idea about relative magnitudes of first- (Oi) and second-order (O2) perturbative corrections appearing in 
our series expansions, we quote their ratio O1/O2 (at fixed point) in the second column. The larger (absolute) values 
of this ratio correspond to the better apparent convergence of truncated series. Except for Pi and 711, the values of 
Oi I O2 are positive. This means that the sighns of the first- and second-order corrections do not alter for most of the 
exponents. Note, that a very similar situation has been encountered in the analysis of perturbation expansions for 
surface exponents at the ordinary transition in pure systems [pj] . 



6 



In fact, as it was shown in [ p4| , the best numerical estimates for surface exponents in pure semi-infiinite systems 
have been given by diagonal [1/1] Pade approximants. Since in the present case the qualitative behavior of underlying 
series expansions is very similar, we expect to obtain the numerical results of comparable rather good quality from 
nearly-diagonal two- variable rational approximants of the types 

[11/1] = l + ^i^ + 5i. + «n^^ (5,5) 

and 

1 -f biu + Oiv + Oiiuv 

The corresponding values of surface critical exponents are given in the last two columns of the Table I. As we can see, 
these numbers do not differ significantly between themselves. 

The values [0/0], [1/0], and [2/0] are simply the direct partial sums up to the zcroth, first, and second orders, 
respectively. We consider the [11/1] and [1/11] values as the best we could achieve from the available knowledge 
about the series expansions in the frames of the present approximation scheme. Their deviations from the other 
second-order estimates of the table might serve as a rough measure of the achieved numerical accuracy. 



VI. 4 - e EXPANSION 



An alternative approach to calculate the desired surface critical exponents is the treatment of the theory in 4 — e 
space dimensions {^/e expansion) and subsequent extrapolation to e = 1. This approach was initiated by Ohno 
and Okabe [|35[| in 1992. These authors considered the two-loop approximation for correlation functions in random 
semi-infinite space. They derived the corresponding series expansions in powers of renormalized coupling constants 
and £, for surface exponents r]\\ and r]±. We quote their results for the ordinary transition in the case of present 
interest, setting n = 1 in the expressions of Ref. [Q and changing normalizations of coupling constants {u — > u/24, 
w — > — w/3) to fit with our notations: 

,||=2---- + - + -.^-H-u. + ... (6.1) 

w w 5 , 11 9 5 „, 
^- = 1-6-4 + 36^ +48^ +T2"" + ---- ^'''^ 

In the present problem dots represent less important terms of order 0{£^). Unfortunately, only the first non-trivial 
corrections oc-^/e have been obtained from these equations in ]3^ . 

Actually, it is possible to derive one more term in the ^/e expansion of surface critical exponents using the known 
expressions for the fixed-point values up to 0{e) p^ , p^ 

/fe , ^^ 110 + 63C(3) ^ 
" =-^V53+'^ 53^ ^^-^^ 



where ({3) ~ 1.2020569 is the Riemann ^-function, and the usual geometric factor Kd = 2^ '^tt '^^'^ /r(d/2) has been 
absorbed into the redefinitions of the coupling constants. Thus we obtain 

[6e 756C(3) - 5 

.. = l~i#+^^ffl^.. (6.6) 

' 2 V 53 2 • 532 ^ ' 

It can be easily verified that the above exponents satisfy the scaling relation = (f? + f?||)/2 with the correct value 

^ = -l^ + 0(ei) (6.7) 
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of the bulk theory. 

Taking into account scahng relations for surface critical exponents and \fe expansions for random bulk exponents 
V and J] [p^p3| we obtain, in addition, 



1 1 /fe 3 252C(3) - 37 

2 ^ iV 53 ^ 8 53 
1 1 /6e 6211-1512C(3) 



711 = -o + 7\/ ^ - o (6-8) 




2 V 53 ' 8-532 

7 (108C(3) - 137) 
16 53" 



e, (6.9) 
2 £. (6.10) 



Similarly as in the case d = 3 of the previous section, we performed a Pade analysis of our expansions at e = 1. 
The numerical values of surface critical exponents obtained in this way are represented in the Table II. 

The [1/0] values for the exponents r\\\ and 77 reproduce the first-order results by Ohno and Okabe On the 

other hand, the other exponents, 711, and 71 slightly differ. The reason is that we calculated our [1/0] estimates 
directly from each respective \fe expansion, while in Ref. ]|35[] they were obtained from the scaling relations using the 
above numerical vales of r\\^ and 77^. 

Comparing the results from Tables I and II we see that the values of first-order approximants, denoted by [1/0] and 
[0/1] in both cases, are of comparable magnitudes. But, on the other hand, the values from second-order approximants 
are significantly different in both tables. The reason is that the second-order contributions of the \fl expansion provide 
corrections of opposite sighns, as compared to the "three-dimensional" theory of the previous section (corresponding 
orders' ratios O1/O2 have opposite sighns). 

Our choice will be in favor of our estimates derived directly in three dimensions for the following reasons. It is well 
known that to the order of 0(e) the _\/e expansion fails to yield the positive correlation function exponent 77 for the 



random bulk Ising system [see Eq. (6.7)]. In fact, if we try, using the scaling relation 77 = 2r7_L — r7||, to reproduce the 
numerical value of 77 from our second-order data of Table II, we also always obtain negative values. This deficiency 
is not present in our calculations directly at d = 3. Moreover, there are several reports in the literature on "bad" 
behavior of the ^ft expansion at larger orders, and for other bulk exponents ] |l^ , ^ , p2t . 

At the same time, the fixed-dimension massive field theory appears to give quite regular and reliable results for 



random bulk systems in three dimensions [p0|-p3[, even at rather low orders of perturbation theory |15 5q , p7| , pi7| . 
From our experience, the massive field theory works also well in description of surface critical behavior in pure 
three-dimensional systems [^,0 . 



VII. SUMMARY 



In the present work we studied the surface critical behavior of three-dimensional quenched random semi-infinite 
Ising systems with free plane boundaries. We have calculated the corresponding surface critical exponents of the 
ordinary transition employing two alternative possibilities: (a) using the massive field-theoretic approach directly in 
d — Z dimensions, and (b) performing the \ft expansion about the upper critical dimension d = 4 with subsequent 
extrapolation to e = 1. In this latter calculation we extend, to the order of 0((-ye)2), the previous first-order results 
by Ohno and Okabe [||. 

We performed a rational approximants' (Pade) analysis of the resulting perturbation series expansions in both 
cases attempting to find out the best numerical values of surface exponents in three dimensions. However, the typical 
behavior of perturbative expansions in both calculational schemes appeared to be qualitatively different. 

In the last section we gave some arguments in favor of results obtained in the framework (a): directly at d = 3. 
Thus, the summary for our final numerical values of the surface critical exponents at the ordinary transition in the 
presence of randomness is 

7711 = 1.36, Ai = 0.53, r;_L = 0.74, /3i = 0.88, 
711 =-0.24, 71=0.83, (5i = 2.1, (Jn = 0.69 . 

The above values stem from the last two columns of the Table I. The estimate of 711, for which no approximants 
[11/1] and [1/11] exist, has been derived from the scahng relation 711 = v(\ — 77||), where we used v = 0.68 ]|l5| and 
the above value of r7|| = 1.36. We believe, these values are the best one can derive from all above calculations. 

The values of our exponents characterising the surface critical behavior of semi-infinite quenched random Ising-like 
systems are apparently different from their counterparts of pure Ising systems 7711 = 1.528, Ai = 0.464, 

77^ = 0.779, /3i = 0.796, 711 = -0.333, 71 = 0.769, Si = 1.966, Sn = 0.582. 
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We quantitatively confirm a general expectation that the change of the bulk universality class of a system should 
affect its boundary critical behavior. So, the semi-infinite systems with quenched bulk disorder are characterized by 
a new set of the surface critical exponents. 

We suggest that the obtained results could stimulate further experimental work as well as numerical investigations 
of the boundary critical behavior of disordered systems. 
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TABLE I. Surface critical exponents of the ordinary transition for d = 3 up to two-loop order at the random-fixed point 
u* = -0.60509, V* = 2.39631. 



q 



exp 


o. 


[0/0] 


[1/0] 




[0/1] 


[2/0] 


[0/2] 


[11/1] 


[1/11] 


m 

Ai 

Vi. 
/3i 
711 
71 

5i 
Su 


2.10 
2.35 
2.62 
-3.11 
0.00 
3.71 
1.90 
1.71 


2.00 
0.25 
1.00 
0.75 
-0.50 
0.50 
1.67 
0.33 


1.676 
0.412 
0.838 
0.912 
-0.500 
0.743 
1.847 
0.477 




1.721 
0.443 
0.861 
0.943 
-0.500 
0.821 
1.868 
0.501 


1.522 
0.481 
0.776 
0.860 
-0.380 
0.808 
1.941 
0.561 


1.581 
0.507 
0.800 
0.841 
-0.364 
0.832 
1.968 
0.595 


1.364 
0.530 
0.736 
0.875 

0.829 
2.056 
0.695 


1.358 
0.531 
0.735 
0.876 

0.832 
2.054 
0.691 




TABLE II. Surface critical exponents 


of the ordinary transition from the v'eexpansion. 




exp 




[0/0] 




[1/0] 




[0/1] 


[2/0] 


[0/2] 


[1/1] 


m 
pi 

711 
71 


-2.09 
-2.22 
37.60 
-2.37 
-4.17 


2.00 
1.00 
0.75 

-0.50 
0.50 




1.664 
0.832 
0.792 
-0.416 
0.668 




1.712 
0.856 
0.794 

-0.408 
0.702 


1.824 
0.907 
0.793 
-0.451 
0.628 


1.792 
0.892 
0.793 
-0.457 
0.611 


1.772 
0.884 
0.793 

-0.441 
0.636 
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